Coulomb energy and gluon distribution in the presence of static sources. 
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We compute the energy of the ground state and a low lying excitation of the gluonic field in the 
presence of static quark -anti-quark (qq) sources. We show that for separation between the sources 
less then a few fm the gluonic ground state of the static qq system can be well described in terms of a 
mean field wave functional with the excited states corresponding to a single quasi-particle excitation 
of the gluon field. We also discuss the role of many particle excitations relevant for large separation 
between sources. 
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I. INTRODUCTION 

Recent lattice simulations lead to many new theoretical 
insights into the dynamics of low-energy gluon modes Q, 
0130100 In the quenched approximation aspects of 
confinement emerge from studies of the gluonic spectrum 
produced by static color sources. In the following we will 
focus on the pure gluon dynamics (the role of dynamical 
quarks in the screening of confining gluonic strings has 
recently been studied in 0). 

Lattice studies indicate that with relative separations 
between two color sources, R > 1.7 fm , the ground state 
energy obeys Casimir scaling 0, 0|. This means that 
the spectrum of gluon modes generated by static color 
sources depends on the dimension of the color represen- 
tation of the sources rather than on the N-ality of the 
representation (which is related to the transformation 
property of a representation with respect to the group 
center) |ldj . For example, for two sources in the fun- 
damental representation, lattice computations show, as 
expected, that energy grows linearly with the separa- 
tion between the sources. However, also for sources in 
the adjoint representation (with N-ality of zero), lattice 
produces a linearly rising potential, even though for van- 
ishing N-ality screening is expected to saturate the po- 
tential. Screening comes from the production of gluon 
pairs (glueballs) which vanishes in the limit of a large 
number of colors. Casimir scaling is thus telling us that 
there is, at least in the energy range relevant for hadronic 
phenomenology, a simple, universal (source independent) 
description of the confining string. 

The lattice spectrum of gluonic modes generated by 
sources in the fundamental representation, i.e., a static 
quark-antiquark (qq) pair, has been extensively studied 
in 0, 0- The ground state energy, which as a func- 
tion of the qq separation is well represented by the Cor- 
nell, " Coulomb+linear" potential and the spectrum of 
excited gluonic modes have been computed. The excited 
gluonic modes lead to excited adiabatic potentials be- 
tween the sources in the sense of the Born-Oppenheimer 
approximation with the quark sources and gluonic field 
corresponding to the slow and fast degrees of freedom, 



respectively 0, |^. The gluonic wave functional of 
these modes can be classified analogously to that of a 
diatomic molecule. The good quantum numbers are: 
A = 0(E), 1(11), 2(A), • ■ ■ which give the total gluon spin 
projection along the qq axis, PC = +l(g),—l(u) which 
correspond to the product of gluon parity and charge 
conjugation, and Y = ±1 which describes parity under 
reflection in a plane containing the qq axis. The ground 
state corresponds to Ap C = £+ . The lattice calculations 
show that the first excited state has the II U symmetry 
(for A 7^ 0, Y = ±1 states are degenerate) and thus has 
PC = -1. 

The lattice spectrum of gluonic excitations is well re- 
produced by the bag model [l3L fl4| . The crucial fea- 
ture of the model that makes this possible is the bound- 
ary condition, which requires the longitudinal component 
of the chromo-electric and transverse components of the 
chromo-magnetic field of the free gluon inside the cavity 
to vanish at the boundary of the bag. This results in 
the TE mode with pseudo- vector, J P ' C = quan- 
tum numbers having the lowest energy, which leads to 
the II U adiabatic potential being the lightest from among 
the excited gluonic states in the qq system. In another 
model, the non-relativistic flux tube model [l5| , the 
PC = +1 quantum numbers of the low- lying gluon mode 
result from associating a negative parity and a positive 
charge conjugation to the lowest order transverse phonon 
(unlike that of a vector field). This also results in the 
IL U quantum numbers for the first excited adiabatic po- 
tential. Finally in a QCD based quasi-particle picture 
the intrinsic quantum numbers of the quasi-gluons are, 
J P ' C = that of a transverse vector field [lti 

the first excited adiabatic potential between qq sources is 
associated with a single quasi-gluon excitation and this 
quasi-gluon interacts via normal two-body forces with 
the sources, then, one expects the quasi-gluon ground 
state wave function to be in an orbital <S*-wave, which, in 
turn, leads to the net PC = +1 and the II S symmetry 
for this state. This is in contradiction with the lattice 
data as noted in [l^. The bag model and the flux tube 
model give the right ordering of the spectrum of low ly- 
ing gluonic excitations, even though they are based on 
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very different microscopic representations of the gluonic 
degrees of freedom. 

There are indications from lattice simulations of vari- 
ous gauge models that the adiabatic potentials approach 
that of the flux tube, or better string-like spectrum for 
qq separations larger then R > 3 fm 0], however, the 
situation for QCD is far less clear Q. In particular, for 
large separations between the sources the splitting be- 
tween nearby string excitations is expected to fall off 
as oc tt/R. The lattice results indicate, however, that 
the spacing between the adiabatic potentials is close to 
constant. At distances R < 0.2 fm the flux tube model 
becomes inadequate while QCD is expected to become 
applicable. For example as R — » 0, the Coulomb poten- 
tial between the quark and the anti-quark in the color 
octet is repulsive, and, indeed, the results of lattice cal- 
culations do seem to have that trend. The bag model 
attempts to combine the perturbative and long range, 
collective dynamics by using a free fled theory inside a 
spherically symmetric bag and deforming the bag to a 
string like shape as the separation between the sources 
increases. A self consistent treatment of bag and gluon 
degrees of freedom is, however, lacking. 

Another model which aims at relating the string-like 
excitations at large qq separations with the QC D g luon 
degrees of freedom is the gluon chain model ja, and 
versions thereof [2(j. The model is based on the as- 
sumption that as the separation between the sources in- 
creases pairs of constituent gluons are created to screen 
the charges in such a way that the Fock space is dom- 
inated by a state with a number of constituent gluons, 
which grows with the qq separation. Recently, support 
for the gluon chain model came from lattice studies of 
the Coulomb energy of the qq pair pll |22| . As shown in 
[23| . at fixed R, Coulomb energy bounds the true exact 
(from Wilson line) energy from above. The Coulomb en- 
ergy is defined as the expectation value of the Coulomb 
potential in a state obtained by adding the qq pair to the 
exact ground state of the vacuum, i.e., without taking 
into account vacuum polarization by the sources. The 
addition of sources changes the vacuum wave functional 
by creating constituent gluons as described by the gluon 
chain model. 

In this paper we discuss the structure of the qq state 
in terms of physical, transverse gluon degrees of free- 
dom. In particular, we focus on the importance of con- 
stituent gluons in describing the excited adiabatic poten- 
tials. For simplicity and to make our arguments clearer, 
we concentrate on excited adiabatic potentials of single, 
Ap^ = symmetry. A description of the complete 

spectrum of excited potentials will be presented in a fol- 
lowing paper. Our main finding here is that a descrip- 
tion based on a single (few) constituent gluon excitation 
is valid up to R ~ few fm, with the gluon chain turning 
in, most likely, at asymptotically large qq separations. 
Consequently, we show how the gluon chain model can 
emerge in the basis of transverse gluon Fock space. 

In Section II we review the Coulomb gauge formulation 



of QCD and introduce the Fock space of quasi-gluons. In 
Section III we review the computation of the ground state 
and the excited £+ potentials. There we also discuss the 
role of multi-particle Fock sectors and a schematic model 
of the gluon chain. A summary and outlook are given in 
Section IV. 



II. COULOMB GAUGE QCD 

In the Coulomb gauge gluons have only physical 
degrees of freedom. For all color components a = 
1, • • • ,Ng — 1 the gauge condition, V • A a (x) = 
0, eliminates the longitudinal degrees of freedom and 
the scalar potential, A a ' a , becomes dependent on the 
transverse components through Gauss's law 24]. The 
canonical momenta, II a (x), satisfy [Ilf (x),A^(y)] = 

-iS ab 6' l rj>(V)5 3 {x-y) where $'(V) = ^-V'VVV 2 ; in 
the Shrodinger representation, the momenta are given by 
II a (x) = — iS/SA a (x). More discussion of the topolog- 
ical properties of the fundamental domain of the gauge 
variables can be found in [25j . The full Yang-Mills (YM) 
Hamiltonian with gluons coupled to static qq sources in 
the fundamental representation is given by, 

H = Hq + H Qg + Hqq, (1) 

where Ho is the YM Hamiltonian containing the kinetic 
term and interactions between transverse gluons. The 
explicit form of the YM Hamiltonian, Hq, can be found in 
|24j . The coupling between qq sources and the transverse 
gluons, Hq q , is explicitly given by, 

H Qg = f dxdyp a Q (x)K[A](x,a;y,b)p b (y), (2) 

where pq = /i t (x)T a /i(x)— (-x.)T* a r](x) is the color den- 
sity of the sources with h and -q representing the static 
quark and anti-quark annihilation operators, respec- 
tively; p — — / a (, ci 7 _1 n & (x) J ■ A c (x) is the gluon charge 
density operator and K is the non-abelian Coulomb ker- 
nel, 

r.-w ,n Q 2 f (1 — A)~ 2 (x, a: z, b) , , 

K A x,a;y,6 - f- / d^- 1 ' ' > , 3 

Air J |z — y| 

with the matrix elements of A given by (1— A)(x, a; y, b) = 
S ab S 3 ( X - y) - s / ac6 V„(l/|x - y|)A c (y). The Faddeev- 
Popov (FP) operator, (1 — A), determines the curvature 
of the gauge manifold specified by the FP determinant, 
J = det(l — A). Finally, the interaction between the 
heavy sources, Hqq, is given by 

Hqq = ^J dxdyp Q (x)K[A}(x, a; y, b)p b Q (y). (4) 

The Coulomb kernel is a complicated function of the 
transverse gluon field. When Hq 9 and Hqq are ex- 
panded in powers of the coupling constant, g, they lead 
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to an infinite series of terms proportional to powers of 
A. The FP determinant also introduces additional in- 
teractions. All these interactions involving gluons in 
the Coulomb potential are responsible for binding con- 
stituent gluons to the quark sources. 



A. Fock space basis 

The problem at hand is to find the spectrum of H for 
a system containing a qq par, 



H\R, N) = E N (R)\R,N). 



(5) 



In the Shrodinger representation, the eigenstates can be 
written as, 

\R,N) = J D[A a (x)}J[A]*%[A a (x)}\*z,i,-Xz,j-A), 

(6) 



with 



|z,j;A) = / l l(|z)r ? ](-|z)|A) (7) 



describing a state containing a quark at position Rz/2 
and color i and an anti-quark at position —Rz/2 and 
color j. We keep quark spin degrees of freedom im- 
plicit since, for static quarks, the Hamiltonian is spin- 
independent. The eigenenergies, En(R), correspond to 
the adiabatic potentials discussed in Section I with N 
labeling consecutive excitations and spin-parity, Ap C , 
quantum numbers of the gluons in the static qq state. 

The vacuum without sources, denoted by |0), in the 
Shrodinger representation is given by, 



|0) = J D[A°(x)]jL4]v|>o[A (x)]|A>, 



(8) 



and satisfies Hq\0) — E vac \0). 

The eigenenergies, E^(R), in Eq. (JSJ contain contribu- 
tions from disconnected diagrams which sum up to the 
energy of the vacuum, E vac . In the following, we will 
focus on the difference, Ejsr(R) — ► £W(i?) — E vac , and 
ignore disconnected contributions in the matrix elements 
of H. 

Instead of using the Shrodinger representation, it is 
convenient to introduce a Fock space for quais-particle- 
like gluons [H EHH H<|. These are defined in the 
standard way, as excitations built from a gaussian (har- 
monic oscillator) ground state. Regardless of the choice 
of parameters of such a gaussian ground state, the set 
of all quasi-particle excitations forms a complete basis. 
We will optimize this basis by minimizing the expecta- 
tion value of the Hamiltonian in such a gaussian ground 
state. We will then use this variational state to repre- 
sent the physical vacuum and use it in place of |0) and 

[A] . The unnormalized variational wave functional is 
given by, W„[A] = (A|0), 



*o[A]=exp (-\ J ^A a (kM|k|)A a (-k) 



(9) 



where A a (k) = J dxexp(— ik-x)A a (x) and the gap func- 
tion, w(|k|) plays the role of the variation parameter. The 
computation of the expectation value of Ho in Wo given 
above, was described in Ref. |2^|. In the following we will 
summarize the main points. 

The expectation value of (0|i?o|0) can be written in 
terms of functional integrals over Z?[A a (x)] with the mea- 
sure J\A\. The functionals to be integrated are products 
of Hq = Hq(TI, A) and the wave functional ^[A]! 2 . For 
example the contribution to (0|i?o|0) from the g = com- 
ponent of the transverse chromo-magnetic field density, 
{B 2 ) = (0| J dx[B a (x)] 2 |0)/(0|0), is given by, 



(B 2 ) = J DAJ[A][V x A a (x)] : 



: *g[A] 
(0|0) 



= N 



dk 



(2tt) 3 2fi(|k| 



(10) 



where Af = 2 x (N 2 — 1) x V counts the total (infinite) 
number of gluon degrees of freedom in volume V and £1 
is the instantaneous gluon-gluon correlation function, 



J DAJ[A]A a (p)A b ( q 



*§[A] 

(0|0) 2fi(|p|) 



(2^ 3 <5(p + q). 

(11) 

In the limit J — > 1, fl becomes equal to the gap func- 
tion u> [261 l2q . Evaluation of functional integrals over 
non-gaussian distributions, like the one in Eq. for 
J 7^ 1 can be performed to the leading order in Nc by 
summing all planar diagrams. This produces a set of cou- 
pled integral (Dyson) equations for functions like Q(p). 
The Dyson equations contain, in general, UV divergen- 
cies. To illustrate how renormalization takes place, let 
us consider expectation value of the inverse of the FP 
operator, 



5 ab d(* - y) = J DAJ[A]g(l - A)" 1 (x, a; y, b) 



*§[A] 



(Q|0) ■ 
(12) 

From translational invariance of the vacuum, it follows 
that the integral depends on x — y and the Dyson equa- 
tion for d becomes simple in momentum space. Defining, 
<i(x — y) — * d(p) = J dxexp(— ik • x)rf(x), one obtains, 
(p = |p|, etc. ), 



1 



1 



dip) 5(A) 



Nc 
2 



dq (1 - q ■ p) 

(27r)an(|p-q|)q: 



d(q). (13) 



As expected from asymptotic freedom, for large mo- 
menta, ft(k)/k — > l + 0(logfc); the integral in Ea.ll2lbe- 
comes divergent as q — > oo, and we need to introduce an 
UV cutoff A. The cut-off dependence can, however, be re- 
moved by renormalizing the coupling constant g — > <?(A). 
The final equation for d(p) , renormalized at a finite scale 
/i, is obtained by subtracting from Eq. I|12ll the same 
equation evaluated at p = //. 

One also finds that the expectation value of (1 — A) 2 , 
which enters in the Coulomb kernel, if [A], requires a 
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multiplicative renormalization. We define the Coulomb 
potential as, 

J DAJ[A]K[A] (x, a; y, = -*«&Vfc(x - y), 

and introduce a function / by, 



V c (k) = J dxe tk x F c (x) = - 



f{k)d\k) 



k 2 



(15) 



This function then satisfies a renormalized Dyson equa- 
tion, 



f(k) 



+ 



m 

Nc 
2 



rig (l-q-p)rf 2 ( g )/( g ) 
(2tt)3 n(|p-q|)q2 



-(* 



M) • 
(16) 



Finally, the bare gap equation, 5[(Q\H \0)/{0\Q)]/6u(k) = 
0, contains a quadratic divergence proportional to ~ A 2 . 
This divergence is eliminated by a single relevant oper- 
ator from the regularized Hamiltonian, the gluon mass 
term, which is proportional to A 2 J dxA a (x) . The renor- 
malized gap equation determines the gap function w(fc), 
and it depends on a single dimensional subtraction con- 
stant, 

The functions described above completely specify the 
variational ground state, and the complete Fock space 
basis can be constructed by applying to this variational 
ground state quasi-particle creation operators, a a '^(k, A), 
defined by, 



A Q (x) 



n a (x) = 



r/k 



(2tt)3 



[a a (k, A)e(k, A) 



+a Q < t (-k,A)e(-k,A)] e jk ' x , 
-a°' t (-k J A)e(-k,A)] e ikx . 



(17) 



Here e represent helicity vectors with A = ±1. This 
Fock space and the corresponding Hamiltonian matrix el- 
ements depend on four parameters (renromalization con- 
stants), u>(fi), d(n), f(/j,) and one constant needed to reg- 
ulate the FP determinant. The FP determinant enters 
into the Dyson equation for O(fe). 

In principle, if the entire Fock space is used in building 
the Hamiltonian matrix and no approximations are made 
in diagonalization, the physical spectrum will depend on 
the single parameter of the theory i.e the renormalized 
coupling (or d(fx), c/Eq. JED- In practical calculations, 
the Fock space is truncated and this may introduce other 
renormalization constants. Goodness of a particular ba- 
sis, for example the one built on the state given in Eq. JjJJl, 
can be assessed by studying sensitivity of physical observ- 
ables to these residual parameters. 



For example, if we define the running coupling as, 
a(k) ee f(k)d 2 {k), so that V c (k) = we will find 

that for large k, a(Jfc) oc (1/ log c (fc))[l + 0(1/ log(fc))] 
where c ~ 1.5, 29], while in full QCD the leading log has 
power c = 1. The discrepancy arises because we used the 
single Fock state, |0) in definition of Vc (and a). This 
omits, for example, the contribution from the two-gluon 
Fock state, as shown in Fig.^ This two gluon intermedi- 
ate state clearly impacts the short range behavior of the 
Coulomb interaction, but, as discussed in |29j |. it is not 
expected to affect the long range part (partially because 
the low momentum gluons develop a large constituent 
mass) . Similarly, in [23 , the role of the FP determinant 
has been analyzed, and it was shown that it does not 
make a quantitative difference leading to Sl(p) ~ 

This is in contrast, however, to the results of |27| . 
We think this discrepancy originates from the difference 
in the boundary conditions which in j27| lead to f(k) = 
1. This makes possible for the gap equation to have a 
solution for uj(k) which rises at low momenta. If f(k) ^ 
1 and, in particular, if f(k) grows as k — ► 0, which is 
necessary if Vc (R) is to grow linearly for large R, we find 
that u>(k) has to be finite as k — > 0. A more quantitative 
comparison is currently being pursued. We also note that 
lattice simulations [3(j are consistent with the results of 
ii. 



In the following, we will thus set J = 1, which makes 
Q = u), and use the solutions for f{k), d(k) and ui{k) 
found in Ref. [29|. 

Finally, we want to stress that the Coulomb potential, 
defined in Eqs. (|14fl . (|15fl . gives the energy expectation 
value in the state obtained by adding the qq pair to the 
vacuum of Eqs. iJSJ, (J5J), i.e, 



(qq\H\qq) = C F V C (R) - C F V C (0), 



with CfVc(0) originating from self-energies, and 



(18) 



\qq) = |i?,JV = 0,FJ+) 
= J— fe f I — z I //' 



R 



10) 

(o|o>- 



(19) 



The state \R,N = 0, £+) refers the the ground state 
(N = 0) with spin-partiy quantum numbers Ap C = 
0(S)+. The energy CfVc(R) should be distinguished 
from Eq(R) in Eq. (J5J. The latter is evaluated using 
the true ground state of the qq system while the former 
is evaluated in a state obtained by simply adding a qq 
pair to the vacuum. Since a qq pair is expected to polar- 
ize the gluon distribution , these two states are different. 
Furthermore, in this work, the \qq) state in Eq. I|19l) is 
obtained by adding the qq pair to the variational state 
of the vacuum and not to the true vacuum state in the 
absence of sources. 
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FIG. 1: The 0(g 4 ), one loop diagrams contributing to the 
leading log term in the expansion of the /3-function in YM 
theory with heavy sources, a) anti-screening dressing of the 
Coulomb potential by transverse gluons, b) Debye screening 
of the Coulomb potential by transverse gluons. The Coulomb 
potential is represented by the dashed line, and sources by 
thick lines. 




FIG. 2: Comparison between V(R) = C F V c {r) from Eq. J2§ 
(solid line) and V(R) = Eo(R) lattice data from (ro = 
1/450 MeV^ 1 ). 



B. Fitting the Coulomb Potential 



As discussed above, the Coulomb energy, CfVc(R), 
represents the expectation value of the Hamiltonian in a 
particular qq state (given in Eq. I|19|) ). which is not the 
same as the true eigenstatc of the Hamiltonian for the 
qq system as defined in Eq. (JjJJ. The latter has energy 
E (R). 

According to 1281 . CfVc(R) > E (R) and numer- 
ical results in further indicate that for large R, 
C F V C (R) ~ a c R and E (R) ~ oR with the Coulomb 
string tension, ac, being approximately three times 
larger then a. In j^j we, however, fitted d(n), /(/i) 
and so that CfVc(R) — > E (R), and a number 
of phenomenological studies have been successful with 
those parameters 0, 1^3, EM 1^3 . It should be noted, 
however, that the results from [22] for CpVc{R) may 
not directly apply to our analysis since the qq state used 
here to define Vc(R) may be different from the one used 
in lattice computations of Vc{R)- Guided by the suc- 
cesses of the phenomenological applications of our ap- 
proach we proceed with fitting CfVc(R) to Eq(R). It 
is clear, however, that since the qq state of Eq. 119(1 
is a variational state, CfV(R) should be greater than 
Eo(R) ,23]. We will nevertheless proceed with the ap- 
proximation CfVc{R) = Eq(R) and examine the conse- 
quences afterwards. 

In |2^], we have found that the numerical solutions 
to the set of coupled Dyson equations for d(k), f(k) and 
uo{k) can be well represented by, 



d(k) 



3.5 



3-6 (TST) 

log(2.41) 
log(fcVm^+1.41) 



for k < 
for k > m 



(20) 



f(k) 



1.41 ft) 



0.97 



1 41 ( '°g(l-82) 

^1(^/^+0.82) 



0.62 



for k < 
for k > 1 



'<)■ 



to(k) = I mg k<m 9 

v ' ' k for k > m a . 



(21) 



(22) 



The parameter m g — 650 MeV effectively represents the 
constituent gluon mass. It should be noticed, however, 
that u>(k) is the gap function and not the single quasi- 
particle energy. This energy, denoted by E g {k) is given 
by, 



Nc 
2 



W7 C 



Vfc(k-q) 



1 +kq 

2w(?) 



(23) 



Since Vc(k) = —f(k)d 2 (k)/k 2 , which for small k grows 
faster then fc 3 , the integral in j Eq. ffify is divergent. This 
IR divergence is a manifestation of the long range nature 
of the confining Coulomb potential which removes sin- 
gle, colored excitations from the spectrum. As will be 
explicit in the examples studied later, residual interac- 
tions between colored constituents in color neutral states 
cancel such divergencies and result in a finite spectrum 
for color neutral states. In the following analysis, we 
will also need the Coulomb potential in coordinate space. 
We find it practical to approximate the numerical FT of 
Vb(k-q) by, 



Vc(r) = br 



rlog c [(rA)- 1 + a] 



(24) 
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with b = 0.20 GeV 2 , a = 0.83, A = 0.63 GeV, a = 1.24 
and c = 1.51. Comparison between CpVc(R) and E (R) 
obtained from lattice computations is shown in Fig. 

We now proceed to the main subject of this paper, 
namely to investigate the difference between Eq (R) com- 
puted using the single Fock space approximation to the qq 
state (i.e without modification of the gluon distribution) 
and the solution of Eq. JSJ which accounts for modifica- 
tions in the gluon distribution in the vacuum in presence 
of qq sources. We will also compute the first excited po- 
tential with the £+ symmetry. 



III. ADIABATIC POTENTIALS 

To diagonalize the full Hamiltonian in the Fock space 
described above, in principle, requires an infinite number 
of states. In the zeroth-order approximation, Eq(R) — 
CfVc{R), a single state with no j quasi-gluons was used. 
At vanishing qq separation, we expect the wave function 
of the system to be identical to that of the vacuum, and 
the approximation becomes exact. One also expects that 
the average number of quasi-gluon excitations in the full 
wave functional of Eq. © increases with the qq sepa- 
ration. We thus start by examining the approximation 
based on adding a single quasi-gluon and truncate the 
Hamiltonian matrix to a space containing \qq) and \qq~g) 
states, 



(qq\Hqq)^ (qq\H\qqg) 
(qqg\H\qq) (qqg\H\qqg) 



\mg) 



E N (R) 



\iig) 



(25) 



The \qq) state is given in Eq. (|19|) . In the quasi-particle 
representation the state with a single gluon and Ap C 
quantum numbers, \qqg) — \R, n, Ap C ) is given by, 



\R,N,Al c ) 



E 



2j + l 



elk 



8wC F N c J (2tt) 3 



£>i*;(k) + VYD%(k)] ^(k)xU\R,K A), 



(26) 



for A ^ and, 

\R,N,Q? C ) 



E 



2.7 + 1 



inCpNc 



■D5'"(k)^(fc) X * A |fl,k,A), 



for A = (S potentials) where 



(27) 



(28) 



a) 



R/2 



b) 



( -R/2 

c) 



nnmhsisinr^ 
e) 

< 



< 

d) 

nsinnrsirsir^ 
* 



f) 



FIG. 3: Matrix elements, (R, k', X'\H\R, k, A). Diagrams a) 
and b) represent gluon and quark self energies, respectively. 
Diagrams c) and d) represent the Coulomb interaction, Vc 
between the gluon and one of the quarks and between the 
two quarks, respectively. In the bottom row, diagrams e) and 
f) describe matrix elements of the interaction term resulting 
from expansion of the Coulomb kernel K[A] in up to one 
power in gluon field. 



and a f = a a ^T a . In Eqs. HU,©, j g is the total an- 
gular momentum of the quasi-gluon. For vanishing sep- 
aration between the quarks, the system has full rota- 
tional symmetry, and j g becomes a good quantum num- 
ber. In general, the system is invariant only under rota- 
tions around the qq axis. It is only the projection of the 
total angular momentum, A, that is conserved and states 
with different j g become mixed. The wave function x^ x 
represents the two possibilities for the spin-oribt coupling 
of given parity, (j g = L g or j g = L g ± 1). It is given by 
Sn\/V2 for £ = 1 and XS^x/y 7 ! for £ = —1, correspond- 
ing to TM (natural parity) and TE (unnatural parity) 
gluons, respectively. Finally rjy determines the behavior 
under reflections in the plane containing the qq axis, i.e., 
the F-parity. 

The radial wave functions, tp 3 ^ (k), labeled by the radial 
quantum number N and j g , are obtained by diagonalizing 
the full Hamiltonian in the Fock space spanned by the qqg 
states alone, i.e by solving the equation, 



PHP\R,N,Al c ) 



VSZ(R)\R,N,M C ). 



(29) 



Here P projects on the \qq~g) state and V^ q ^,{R) are the 
bare energies of the excited adiabatic potentials, i.e., 
without mixing between states with a different number of 
quasi-gluons. Analogously, CpVc(R) is the bare ground 
state energy Eq(R). The matrix elements of PHP are 
shown in Fig. [3] and given explicitly in the Appendix. 
The mixing matrix element, 



(qq\H\qqg) = v%g 9s (R), 



(30) 



depends on the number of bare, qqg states from Eq. I|29|) 
kept, N = ■ N max and the separation between the 
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# 



FIG. 4: The matrix element (R\H \R, k, A). It originates from 
expansion of the Coulomb kernel K[A] to first order in A 

sources, R. It is shown in Fig. 0] and given in the Ap- 
pendix. 

The (N max + l) x (N max +1) Hamiltonian matrix shown 
in Eq. (|2*5|l is explicitly given by, 




C F [V c {R) - V G (0)] N = M = 

Vg™ g9 (R) N = 0,M=1-N max 

N {R) N = 1 - N max , M = 

V™ (R)S NM N,M = 1 — N max 



(31) 

IV. NUMERICAL RESULTS 

In terms of £ and 77, the PC and Y quantum numbers 
of the gluonic field are given by, 

PC = e(-D-+\ Y = { ^<"£ A A f °; A Q * ° . (32) 

In the following, we will concentrate on the states with 
A = 0, PC = g{+) and Y — +, i.e., of £+ symmetry, 
since it is only these states that mix the bare \qq) state 
with the states with non-vanishing number of gluons. 

For the £+ potentials, the wave function contains TM 
gluons, £ = 1 of natural parity and PC = +1 which im- 
plies j g = 1,3, • • • . As discussed above, for R — > 0, j g 
becomes a good quantum number, and we have verified 
numerically that for R in the range considered here the 
contributions from j g — 3 and higher are at a level of a 
few percent. Diagonalization of the Hamiltonian in the 
qqg subspace alone, leads to the V^ q ^ (R) potential which 
is shown in Fig. [S] (upper solid line) for the lowest exci- 
tation with N = 1. The dashed line is the result of using 
the one- and two-body interactions depicted in Figs.|3Ja,-d. 
(i^jj — ijj^ in Eq. (|39l)). These are also the interactions 
that were used in 17]. When the three-body interac- 
tions shown in Fig. E]3,f are added, the energy moves up. 
This discrepancy is then also a measure of how far our 
variational, truncated Fock space expansion is from the 
true excited state. The three-body potential is expected 




FIG. 5: Comparison between V(R) = CVVc(r) from Eq. J2§ 
(solid line) and the V(R) = Eo(R) lattice data from Q (ro = 
1/450 MeV^ 1 ). 



to be responsible for reversing the ordering between the 
Tl u and II g surfaces; with only one- and two-body in- 
teractions, the Tig potential has lower energy than II U , 
which is inconsistent with the lattice data [ijj]- In the 
Appendix, we also show that the three-body term is sup- 
pressed at large separations, and thus the net potential 
approaches the Casimir scaling Cf^R limit as R — > 00. 
Finally, we note that when the Fock space is restricted to 
single quasi-gluon excitations, the diagrams in Fig.[3]and 
Fig. [35] represent the complete set of Hamiltonian matrix 
elements . 

The general features of higher excitations, V r ( S 9 J ^(i?) 
for A > 1, follow from the structure of the Hamilto- 
nian, which represents a one-body Schrodinger equation 
for the single quasi-gluon wave function in momentum 
space. The kinetic energy corresponds to the one-body 
diagram in Fig. and the potential to the diagrams 
in Fig. 0t,e,f. The diagrams in Figs. Eh>,d give an in- 
dependent shift describing the qq self-interactions and qq 
octet potential. The IR singularity in the gluon kinetic 
energy, E g , is canceled by the collinear singularity of the 
two-body potential, the qq self energy and qq octet poten- 
tial. On average, gluon kinetic energy contributes an ef- 
fective quasi-gluon mass of the order of m g . Quasi-gluon 
are thus heavy, and adding Fock space components, with 
more gluons, \qq, 2g), ■ ■ ■ \qqn g g), for small- R will result 
in higher adiabatic potentials with (N = 2, 3 • ■ • ) that 
are split from the first excited state by ~ n g m g . At 
large R 7 the two-body Coulomb potential dominates and 
together with Coulomb energies of the pair-wise gluon 
interactions, results in the Casimir scaling (we will dis- 
cuss this in more detail in the following section). In the 
absence of mixing between Fock space components the 
number of quasi-particle gluons in the \qq,n g g) state is 
conserved, and they directly map in to the tower of ex- 
cited adiabatic potentials. 
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Lattice 

qq with qqg mixing 
qq without qqg mixing 
qqg, 1st excited 



r/(2r ) 




r/(2r.) 



FIG. 6: Comparison between V(R) = C F V c (r) from Eq. J2§ 
(solid line) and the V(R) = Eo(R) lattice data from Q (ro = 
1/450 MeV" 1 ). 



FIG. 7: Normalized probability of finding the bare \qq) state 
in the full ground state of the \qq, N — 0) (which is also equal 
to the probability of finding the \qq~g) state in the first excited 
\qq,N = 1) state). 



We will now address the effects of mixing between \qq) 
and | qqg) states. The only non- vanishing diagram is 
shown in Fig. [23 Since, as discussed above, the V^, q ^(R) 
potentials are split from the first excited state, N = 1, 
by at least m g , the mixing matrix in Eq. (|31|l saturates 
quickly, and in practice, only the N = 1 state is relevant. 
However, even this single state mixing leads to a very 
small energy shift. In Fig.[S]the dashed line corresponds 
to the energy of the ground state without mixing, (the 
same as the solid line in Fig.EJ), and the solid line shows 
the effect of mixing. The effect of the mixing is small. 
Numerically, we find that the full ground state, 

\qq, N = 0) = Z°(R)\qq) + Z° q \qqg), (33) 



i 

nnnnnn^ 




is still dominated by the \qq) component and the first 
excited state, 

\qq, N = 1) = Z\ q {R)\qq) + Z^qqg), (34) 

by the \qqg) component. The probabilities of each are 
shown in Fig. We see that, for distances between 
sources as large as 5 fen, the admixture of the gluon com- 
ponent is only of the order of 10%. 

This small admixture of the \qq,g) in the full ground 
state is correlated with the small shift in the £+ sur- 
face shown in Fig. and would justify using the ground 
state, exact energy to constrain the Coulomb poten- 
tial Vg. This is, however, contradicting the results of 
Ref. [21j where the effect of mixing must be large since it 
results in a factor of three in the ratio of the unmixed to 
mixed string tensions. One possible explanation is that 
there is an accidental suppression of the mixing interac- 
tion matrix element for the two states considered here, 
|gg) and \qqg). Inspecting Eq. Q51f) . we note that due 
to the gradient coupling of the transverse gluon to the 
Coulomb line, the coupling vanishes both for small and 



FIG. 8: Matrix elements (qq\H\qq, 2g) leading to the \qq, 2g) 
component in the ground state E^j~ potential, a) interaction 
mediated via the Coulomb line coupled to quark sources, b) 
interaction between a single quark and the gluon charge den- 
sity. 



large-R. In contrast, a two gluon state can be coupled to 
|gg) with either the Coulomb line mediated interaction 
as shown in Fig. |HK or the quark density- gluon den- 
sity interaction shown in Fig. As discussed in the 
Appendix, at large distances the former is suppressed 
and it is easy to show that the latter is proportional to 
CpV C (x - R/2) + C F Vc(x + R/2) (once the gluon spin 
is neglected) and persists at large distances. In the large- 
N c limit C F = N c /2(1 + 0(l/N c ))- It is therefore pos- 
sible that the \qq,2g) component of the full |gg, N = 0) 
state is actually more important then the \qq~g) one. We 
will investigate this further in section llV Al 
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R/2 



m 
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FIG. 9: Typical diagrams contribution to mixing between n 
and m gluon states. Vertical dots represent any number of 
gluons not affected by the interaction, a) mixing mediated by 
the Coulomb potential, b) same as in a) with rearrangement 
of gluons, c) long-range Coulomb interaction between gluon 
charge densities, d) same as in c) but with the charge density 
of the quark sources. 



A. Multi-gluons states and the chain model 

As shown above, the quasi-gluon degrees of freedom 
defined in terms of a variational quasi-particlue vacuum 
provide an attractive basis for describing gluon excita- 
tions. This is in the sense that for source separations 
relevant for phenomenology the color singlet states can 
be effectively classified in terms of the number of quasi- 
gluons. This basis, however, does overestimate the ener- 
gies (as expected in a variational approach) , and this fact 
together with lessons from other models can give us guid- 



ance for how to improve on the variational state of the qq 
system. As the separation between quarks increases one 
expects the average number of gluons in the energy eigen- 
state to increase. This is because it becomes energetically 
favorable to add a constituent gluon which effectively 
screens the qq charge. Furthermore, the spacial distribu- 
tion of these gluons is expected to be concentrated near 
the qq axis in order for the energy distribution to be that 
of a flux tube, as measured by the lattice. An improve- 
ment in the ansatz wave functional will therefore result 
in a more complicated Fock space decomposition with a 
large number of quasi-gluons present, even at relatively 
small separations between the sources. In this section we 
will first discuss how multi-gluon states indeed become 
important, even in the case of the quasi-gluon basis used 
here. We then compare with expectations from other 
models and discuss the possible directions for improving 
the quasi-gluon basis. 

As discussed in the Appendix, at large separations the 
interactions between multi-gluon Fock states mediated 
by the Coulomb potential, shown in Fig.|H^,b, require all 
but two gluons to be at relative separations smaller than 
R. Furthermore, rearrangement of gluons leads to 1/Nc 
suppression. For large R, the largest diagonal matrix 
elements of H are the ones corresponding to the long- 
range Coulomb interaction between charge densities as 
shown in Fig. |^;,d. To leading order in Nc, the gluons 
should be path-order along the qq axis. For simplicity, we 
will neglect the gluon spin and use a single wave function 
to represent a state with an arbitrary number of gluons. 
We write 



\QQ,n g g) = N n 



R/2 px„ g « 3 rx 2 

dx n at ! (x n ) dx n lO^ixn i) • • • / dx 2 a^(x 2 ) dxia*(xi)\0) (35) 



R/2 



R/2 



r 



R/2 



R/2 



where we have also forced all gluons to be on the 
qq axis. The factor N rlg = (n g \/C F N c R) 1/2 is, 
to leading order in Nc fixed by the normalization 
condition, {qq,n g g\qq,n' g g} — S ng , n ' g , where we used 

[a(xi), ar(xj)] — Sij. In this basis, the diagonal matrix 
elements of the Hamiltonian (cf. Fig.|5J:,d) add up to 



{qq,n g g\H\qq,n'g) = C F V C (R) C F bRS n 



(36) 

The off-diagonal matrix elements are dominated by in- 
teractions between color charges, e.g., similar to the ones 
in Fig. [HJd, but with the upper vertex attached to a gluon 
line. With the approximations leading to Eq. (|33|) a ver- 
tex which either annihilates or creates two gluons results 
in a vanishing matrix element since in our basis no two 
gluons are at the same point. Smearing each gluon in the 
coordinate space by a distance of the order of 1 /m g will 



give a finite matrix element , which just like the diagonal 
matrix elements grows linearly with R, 



(qq,n g g\H\qq,n'g) 



jC F bR 



Sn„,n'+2 + Sn g ,n'+2 



(37) 



where 7 is a parameter representing the effect of a smear- 
ing, and we expect I7I < 0(1). In addition, each gluon 
has a kinetic energy of the order of m g , so H nn — > 
H n n + nm g . The model Hamiltonian can be easily diago- 
nalized numerically, and in Fig. 1101 we plot the energy of 
the ground state and of the first excited state as a func- 
tion of R. It is clear that in the absence of accidental 
spin suppression, which, as discussed earlier, takes place 
for the (qq\H\qqg) mixing matrix, the effect of the mixing 
with two and more gluons can produce shifts in the lowest 
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' r/(2r„) 



FIG. 10: Shift in the ground state Eg energy due to cou- 
pling with muti-gluon states of the model Hamiltonian of 
Eqs. 1361 . 13711 . The maximum number of states was taken to 
be n g , m ax = 40. The other parameters are b = 0.21 GeV -2 
and m g — 0.65 GeV. 




FIG. 11: Average number of quasi-gluons in the full eigen- 
state of the model Hamiltonian of Eqs. 13611 . 1371 . 



adiabatic potential and decrease the Coulomb string ten- 
sion by as much as a factor of ~ 2 at R ~ 3ro = 2.6 fm. 
Finally, in Fig. ^2 we plot the average number of glu- 
ons in the ground state of the model Hamiltonian. As 
expected, the number of gluons grows with R; however, 
still a small number of quasi-gluons contributes to the 
ground state at these separations, which again provides 
justification for the quasi-gluon description. 



V. SUMMARY AND OUTLOOK 

We computed the ground state energy and the energy 
of the first excited qq potential with the Ap C = sym- 



metry. We used the quasi-particle basis of constituent 
gluons based on a variational ground state to build the 
Fock space representaion. We found that the qq state 
can be well approximated by a superposition of the bare 
qq state and a few quasi-gluons. The exact computation 
in which the bare qq sate mixes with a state containing 
a single quasi-gluon leads to negligible change in the en- 
ergy of the bare (Coulomb) qq system. We found that this 
is due to an accidental small mixing matrix element of 
the Coulomb gauge Hamiltonian. We have discussed the 
general properties of the mixing matrix between states 
with an arbitrary number of gluons, and using a simple 
approximation, we have found a good agreement with 
the lattice data. The lattice data indicates that there 
is a change in slo pe b etween the Coulomb and the true, 
Wilson potential |2l|]. Based on the representation used 
here, we interpret this in terms of quasi-gluon excitations 
rather than in terms of a flux-tube-like degrees of free- 
dom. We also note that lattice data on splitting between 
several excited qq states does not unambiguously show a 
string- like behavior for separation as large as 2 — 3 fm 0] . 
In fact, the splittings are almost constant, although why 
lattice data has such a behavior is not completely un- 
derstood (including a possible systematic error) |35|]. In 
fact this data is consistent with the quasi-gluon picture 
where each quasi-particle adds kinetic energy of the order 
of the effective gluon mass. The full excitation spectrum 
as well as distribution of energy density is currently being 
investigated. 
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VII. APPENDIX 

Here we list matrix elements of the Hamiltonian in 
the basis spanned by \qq) — \R,N = 0,Ap C ) and 
\qqg) = \R,N ^ 0,Ap C ). The \qq) state exists only 
in the Ap C = configuration. Thus mixing matrix el- 
ements are non- vanishing for \qq~g) with £+ spin-parity 
quantum numbers only. 

For each j g , the wave functions ipN.j {k) are expanded 
in a complete orthonrmal basis of functions 4> m . j g (k) 

^N, jg (k)= »» A ^,(*) ( 38 ) 

m— 1 

with normalization, / f^4>* m >,j' g {k)4> m ,j g (k) = 
8m',m,fij' d g - The expansion coefficients are computed by 
diagonalizing the {m m axj g ,max) x (m max j gimax ) matrix, 
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H, 



obtained by evaluating the diagrams in 



H3 = H 3a + H; 



3b 



H. 



3/- 



(39) 



evaluated in the basis of functions <pm,j g ■ I n numerical 
computations for each j gi we used a momentum grid as 
the basis functions. The numerical results presented were 
for a single j g determined from Eq. ll.'?2H after verifying 
that increasing j g changes the computed spectrum by at 
most a few percent. For arbitrary Ap C the Hamiltonian 
matrix elements are given by, 



H 3a = 



J q Jg 



dkk 2 

W) 1 



c/>* m , tj9 (k)E g (k)<t> m , jg (k), (40) 



— —CpVc{G)5 m i^ m 5j^j g 

/dkk 2 
j2nyV c (k)8 m ., m 5 JiJa) (41) 



with 



V c (k) = - 



d 2 (k)f(k) 
k 2 



(42) 



H 3C = ^ E 



0m' . 



2 ^r-^ 



Vc(x-*) + Vc{x+j) 



ix-(k-q) 



v /(2^ + l)(2 Js + l) 



16tt 



M*0 u>(q) 
oj(q) + V ^( fc ) 



Nc 
2 



E 



A, A' ,£T,(7' ,/i 



(2^/ (2^^'.^^^^-^ 



e -»f '(k-q) , g if .(k-q) 



4>,n,j g (k) 



^(2j' g + l)(2 Jg + l) 



16tt 



^IV'(q)^i 9 A*(k)xLxiA'^:(q)^A(k) + W^(A - -A) 



Mfc) , Mg) 
w(g) V ^( fc ) 



(43) 



and 77^ and £ related to j g and Ap C through Eq. l|3*2^l . 



H 3d - ~2j^Vc(R)(> m >,m5j g ,j g 
1 f dkk 2 

(44) 
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v - f dk dp dc 

3e = ^ J J2nf J2nf (2^ 



dk dp dq mjg (fc) 



dxdydz 



) 3 ^2^) ^Mfc) 
A"(x- 5l,z + y-x,y + 5:) + (R^ -R) 



gix-kgiz-qg-iyp 



v /(2j-+l)(2j g + l) 



(2tt)3 (2tt)3 (2tt)3 v ^(fc) 

yj(2j' g + l)(2 Jg + 1) 



(k+p+2q) + (R ^ _R 



-A) 



e 2 



8tt 



^L'(p)^;;(p)xS'v^ ) o(q)^A'(k)^ ) A(k)xi, A ^; *(q) + Vvv'yiA - -A) 



(45) 



where the sum is over /i, v, A, A', cr, cr' and the kernel is Finally, 
given by 



gixkgiypgizq 



^ (X ' Z ' y) = / (2^)3 (2% ( 2 t)3 ^g^) 

(46) 



and 



^ = q2 Nld(k)d( P )d( q ) 



4 k 2 q 2 p 2 
x [d(fc)/(*0 + d(p)/(p) + d(«)/(«)] (47) 



V /" rfk rf P ^'.j^^m,j a (fc) 
3/ (2^)3 (27r) 3 (27r) 3 ^2^) ^2^) 

x J dxdydz 

^J(2j' g + l)(2 Jg + 1) 



#(x - Z + y _ Xj y _ 5) + (R _> _R) 



ix-k iz-q — iy-p 

e e e 



cStt 



i?iV(p)^;;(p)xS^^,o(q)^%*(k)^,A(k)xi A ^(q) + nm'yiA -> -A) 



\ - /" dk dp dq ^m'jfo) (j> m ,j g (k) 
(2tt)3 (2tt)3 (2tt)3 ^2^) T^f) 

y / (2^ + l)(2 Js + l) 



K(k + q, q, p + q) le^-Ck-P) + (R _> _R) 



8?r 



^V'(p)^;;(p)xL^,o(q)^A*(k)^i,A(k)xi A ^(q) + W^(A - -A) 



(48) 



In the large- Nc limit, g^/Nc ~ 0(1), and since d(fe) oc 5 
and / ~ O(l), all of terms above are 0(1) except Hd 
(which corresponds to a non-planar diagram, see Fig.|3J). 
The products of the three factors, d(pi)/pf, originate 
from the three dressed Coulomb lines in diagrams e and 
/ hi 

Fig. |3| and the three factors of / come from the three 



possibilities to insert the V 2 operator on these three 
lines. The derivative coupling between transverse and 
Coulomb gluons leads to the extra q 2 factor in the nu- 
merator in Eq. I|47|l . In coordinate space this implies that 
A'(x, z, y) is short-ranged in z. Furthermore in each of 
the three terms in Eq. I|47|) there is only one combina- 
tion, d 2 (pi) f (pi) I p 2 , which in momentum, space leads to 
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the confining potential Vc ■ The remaining two are of the 
form d(pi)/pf with d(p) oc which for small mo- 

menta also leads to a short-ranged interaction decreasing 
as for large r. We thus conclude that for the 

three interaction lines connecting the four vertices in the 
" three-body force" of Fig. [3^ only one is long-ranged and 
all others are short-ranged. Along these lines one can 
approximate i^(x, z, y) as 



At large separation R with the wave functions peaking 
at |x| ~ 0, we find that H e grows less rapidly than two- 
body interactions. This is in general true for interactions 
originating from the expansion of K[A] in powers of A 
which couple multiple gluons. This is the basis for the 
approximations discussed in Section. IIV Al 



K(x,z,y) oc <5(z) 



m g Vc{x) m g Vc(y) 



l(m g \y\) a K|x|) c 



(49) 



with < a < 1. Ignoring the gluon spin and all spin- 
orbit couplings we then obtain, 



H: 



3e 



dxdy 



fe'( x ) <Pm{y) 
y/2mg~ yj^m g 



R R , / 

^(x--,y-x,y+-) + (R. 



-R) 



oc / dx^,(x) 



^c(x-f) 
(m s |x+f |)« 



(R 



R) 



4>m (x). The off-diagonal matrix element of the Hamiltonian 

mixing the \qq) and qqg) states, shown in Fig. 01 is given 
(50) by, 



dk dq <t> m ,j g (k) 



dxdz 



, R R . . 

tfi(x--,z-x--)-(R- 



-R) 



v^ 9 



■in 



-°AkA'( k )^,A( k )XvA^ :o (q)=^ J (27r)3 (2?r)3 ^g ==Xi(k + q,q) 



Jf .(k+2q) _ (R ^ _ R ) 



vw 9 



4tt 



-D^ y (k)Di A (k) x {, x Dl:;(q) 



(51) 



distances. Approximately, we find 



#i(x,y) 



dp dq 

(2tt) 3 (2tt) : 



#(p, g)e lxp e lyq 



(52) 



with 



tfi(p,?) 



A^CV^F d(p)d(q) 
2 q p2 q 2 i d (P)f(P) 



d(q)f(q)}. 

(53) 

As expected in the large Nc limit K\ = O(l) and just 
like the three-body kernel described previously, K \ (x, y) 
has mixed behavior for large separations. A term, in 
momentum space, proportional to d 2 f in one of the two 
momentum variables leads to Vc in the corresponding 
position space argument. While for the other momen- 
tum variable it leads to a less singular behavior for large 



*.fr-f.x + !)« 'ff-;W + (R~-H)( M ) 
2 2 (rn g \x+f\)P 



with 1 < (3 < 2. In this limit, ignoring spin dependence, 
one finds 



fL - - / / ,/x/v", ( |x - |x+ (55) 



Thus, similar to the case of i?3 e , we find that at large 
separations the mixing terms grow less rapidly with R as 
compared to two-body interactions. 
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